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Abstract. High degrees of intensity correlation between two independent lasers were observed after prop- 
agation through a rubidium vapor cell in which they generate Electromagnetically Induced Transparency 
(FIT). As the optical field intensities are increased, the correlation changes sign (becoming anti-correlation). 
The experiment was performed in a room temperature rubidium cell, using two diode lasers tuned to the 
*^Rb D2 line (A = 780nm). The cross-correlation spectral function for the pump and probe fields is numer- 
ically obtained by modeling the temporal dynamics of both field phases as diffusing processes. We explored 
the dependence of the atomic response on the atom-field Rabi frequencies, optical detuning and Doppler 
width. The results show that resonant phase-noise to amplitude-noise conversion is at the origin of the 
observed signal and the change in sign for the correlation coefficient can be explained as a consequence of 
the competition between EIT and Raman resonance processes. 

PACS. 32.80.Qk Coherent control of atomic interactions with photons - 42.50.Gy Effects of atomic co- 
herence on propagation, absorption, and amplification of light; electromagnetically induced transparency 
and absorption 

1 Introduction 



^ Permanent address: Departamento de Fisica, Facultad de 
Ciencias Fi'sicas y Matematicas, Universidad de Concepcion, 

Av. Esteban Iturra s/n. Barrio Universitario Concepcion - Electromagnetically Induced Transparency (EIT) has re- 
Chile ceived great attention in recent years in connection to sev- 
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eral interesting phenomena, such as light storage and slow analysis frequency. As the laser intensities are increased, 
light propagation [ll2t3] . The strong interaction between intensity correlations become anti-correlations. Correla- 
light and material media in this situation has been the tions and anti-correlations as big as 0.65 and -0.65, re- 
source of inspiration for various proposals of applications spectively, were observed. The results can be explained in 
of EIT to the quantum manipulation of information and to terms of the conversion of Phase-Noise into Amplitude- 
transfer coherence from light to an atomic medium |2l4l5| . Noise (PN-to-AN) as the lasers interact resonantly with 

rr^i, •, i-ui J -L. n Tj thc atomlc mcdlum. Thc passagB from corrclatiou to autl- 

i he strong interaction between pump and probe nelds ° 

■ -c^T^r -4-1, • 4- correlation can be seen as a consequence of the passage 

m hi 1 can lead to signmcant changes m the noise spectra ^ f b 

r ^ ■ , 1 i 1 ri i- • i • from EIT to a Raman resonance, both present in a 3-level 
oi two independent lasers after propagation m an atomic 

1, . • 1 i- u i 4-1, c ij atom in the A configuration, 

vapor, resulting m correlation between the helds, as was ° 

first observed in Ref. 6 . A recent paper [7| reported an t^,^ 4- a at • • ^ • i u ^ j 

f 1^ ir-i f PN-to-AN conversion m atomic vapors has been stud- 
experimental investigation of the dependence of the phase . 1 . ir>r>i 1 -i T i- i 1 • 1 1 i- 
^ Off since 1991, when its application to high-resolution 

correlation between both fields generating EIT as a func- ^ r. ^ i-j.iirn 

^ spectroscopy was hrst suggested by Yabuzaki et a/.[9j. 

tion of the optical depth and transparency frequency win- mi • • ^■ ^^ ^ ^ ■ ^- ri- i i 

^ ^ f J M. l^jg conversion relies on the characteristics of diode lasers, 

dow. The phase modulation (PM) introduced on the pump i ■ i i -i i ■ r i ■ ^ i 

^ ' -t r- which exhibit excess phase noise tor usual experimental 

field could be read in the probe field signal after inter- h n i i 1 1 4-1 j • u 

° conditions [ll)lll] while the amplitude is generally very 

action with the atoms, for a Raman detuning within the ^ 1 1 mi • ^ ■ ^ i-i i ■ 

' ^ stable. Ihis excess phase noise generates amplitude noise 

transparency window. Nevertheless, no distinction between r ,1 1 • 4_- • j j • 4_i ^ • j- rmrTn 

^ ' of the polarization induced m the atomic medium |9 1 12] 

positive and negative correlation was reported in that pa- j ■ 4_ a 4_ 4_- ^ r ^ 

^ ° t- f Qj^f-i^ ^j-^g intensity nuctuation spectrum ot the transmit- 

per. The possibility of negative correlations for large Ra- ,iti;-4. ij . ^1,1 r •j4.i,rToi 

^ ^ ted light IS strongly dependent on the laser linewidth|13j. 

man detunings is discussed in Ref. 8 , which presents the- o- ^^ n t ^ ^ ^ ^ ^ ^^ t ■ ^ 

- -f Since the nelds detected alter the medium are given by 

oretical calculations for the correlation between the pump ,1 • ■ n u 1 ,1 •, j i ■ i- 4_i j 

'- the input nelds plus the excited polarizations, they end 

and probe fields in EIT configuration caused by the atomic . . . . t-, , 1 , 

^ up acquiring excess noise m the amplitude quadrature. 

dipole fluctuations. mi ^ 1 • ^ ^^ ^ ^ ^ ^ ■ 

Ihe spectral noise components that match atomic res- 
in this paper, we report new measurements in which onance frequencies present larger amplitude oscillations, 
this kind of correlation originates from two independent In this way, it is possible to acquire information about 
lasers. Both fields excite an atomic sample forming a A the medium from the power spectrum of light after the 
system, resulting in a EIT situation. We define a normal- sample [5ll4ll5 j. The influence of laser fluctuations on the 
ized correlation coefficient C, bounded by -1 and -1-1, and atomic polarization has been widely studied for two-level 
report measurements of C as a function of intensity and svstems [16ll7ll8| and many models have been proposed 
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for treating the field phase fluctuation [19|12j . The phase tern as 3-level atoms at rest, agreement with experimen- 
diffusing model has received more attention owing to its tal data is considerably improved by integrating over the 
proximity to the diode lasers extensively used in labora- atoms' different velocity classes and including all the rele- 
tories. In a recent experiment [20], we performed measure- vant excited atomic levels. We also found that the optical 
ments of intensity noise spectra between the a"^ and detuning is essential for explaining the change in sign for 
components of a single, linearly polarized, exciting field in the correlation coefficient. For a perfectly resonant A sys- 
a Rb atomic sample, as a function of the magnetic field tem, the EIT process is dominant and the fields become 
in a Hanle/EIT configuration. We observed correlations as correlated but, for an optical detuning of the order of the 
well as anti-correlations between the different polarization excited-level decay rate, the Raman process prevails and 
components, depending on the detuning, controlled by the the fields become anti-correlated. 
magnetic field. A similar experiment was performed by 

Sautenkov et a/.[2IJ, using two initially phase-correlated 2 Experimental setup 

beams in time domain, who also explained it in terms 

Our experimental setup is shown in Fig. [TJ We employed 

of PN-to-AN conversion [HI . In Ref. [53] PN-to-AN con- 
two external-cavity diode lasers (ECDLs) of i MHz linewidth 

version has also been identified as a source of frequency 

and 15 mW power after optical isolators, tuned to the Ru- 

instabilities in Rb atomic clocks, and was eliminated with 

bidium D2 fine (A — 780 nm). The two beams had linear 

a buffer gas cell that broadens the resonances through col- 
orthogonal polarizations and were combined in a polariz- 

lisions. 

ing beam splitter (PBS). Their powers were adjusted to 
In the present work, the emphasis is on the PN-to-AN have equal intensities at the vapor cell. A small portion 
conversion as a source of correlation between initially in- of Laser 2 was sent to a saturated absorption setup for 
dependent macroscopic fields. The paper is organized as fine tuning, and had its frequency locked to the cross-over 
follows: in section 2 we describe our experimental setup peak between the 5Si/2{F = 3) ^ ^Pi/2{F' = 2) and 
and in section 3 we present a theoretical model that in- 55i/2(^' = 3) ^ 5P3/2{F' = 4) transitions of *^Rb. The 
eludes two phase diffusing fields interacting with an atomic rejected output of the polarizing beam splitter is used to 
system. In section 4 we present our results, beginning with observe EIT in an auxiliary vapor cell. Laser I, tuned to 
the experimental correlation spectra as a function of the the 5Sif2{F = 2) ^ 5^3/2(-P"') transition, is then locked 
analysis frequency and optical intensity. Then, in section on the EIT resonance using a Lock-in amplifier (see Fig. 
4b, we present the numerical results and discussion. As de- [1}. In this way, we guarantee that the Raman resonance 
scribed below, although it is possible to extract the basic condition for EIT was always fulfilled. Both lasers were 
aspects of the phenomena by modeling the atomic sys- locked only by feedback applied to their external cavity 
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After the cell, the beams were separated at a second 
PBS and then analyzed at two independent balanced- 
detection schemes. Photocurrents were combined in ac- 
tive sum/ subtraction circuits (SD), and noise was mea- 
sured with a Spectrum Analyzer. Effective bandwidth of 
our detection is limited by the gain of the amphfiers in 
the range of 2.5 to 14 MHz. Beyond this frequency, elec- 
tronic noise reduced the resolution of our measurements. 
We can, therefore, measure the Sum Ss(lo) and Difference 
SdiLo) noise spectra, as well as the individual laser noise 
spectra 6*11(0;) and 622(0;) by blocking a beam on each 
balanced detection. It is then possible to obtain the nor- 
malized correlation coefficient defined by 

612(0;) 



C(o;) 



(1) 



^6'ii(o;)622(o;) 
The symmetrical cross-correlation spectrum Sij{uj) be- 
tween the lasers i and j is defined by 



Fig. 1. (a) Experimental setup. The saturated absorption 
setup for laser 2 is not shown. 01: optical isolator; PBS: po- 
larizing beam splitter; PD: photo-detector; NF: neutral filter; 
A/2: half-wave plate. (b) Energy level representation for our 
atomic model and its correspondence to the relevant Rb hyper- 
fine states. Solid line represents the pump field with frequency 
oji; dashed line represents the probe field with frequency uj2. 

gratings. The laser intensities were controlled by neutral 
density filters inserted just before the main vapor cell. We 
also used a 2 mm-diameter diaphragm to spatially filter 
the laser beams, ensuring a good spatial superposition and 
a fiat, nearly top-hat, intensity profile over the cell. 



S.jiio) = -{5I,{iu)5I,{u;)* + 5I,{oj)5I,{iu)*). (2) 



Thus, 612 can be obtained from 



Ss — 611 + 622 + 26*12 J 



Sd — 611 + 622 — 26*12 ■ 
6*12 = ^(5*5 — Sd) 



(3a) 
(3b) 
(3c) 



A summary of all possibilities allowed by our setup is 
presented in Table [TJ If Laser 1 is blocked, one uses the 
SD2 circuit to calibrate shot noise (with the subtraction 
position of SD2) or measure the total noise (with the sum 
position of SD2) of Laser 2. If one blocks the Laser 2, an 
analogous reasoning is valid for the SDl circuit. The sum 
and difference noise spectra expressed in ([3|) are obtained 
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Measurement 


Beam blocked 


SDl 


SD2 


SD3 


Total noise 1 (S'li) 


Laser 2 


+ 


n.i. 


n.i. 


s.n. of Laser 1 


Laser 2 




n.i. 


n.i. 


Total noise 2 {S22) 


Laser 1 


n.i. 


+ 


n.i. 


s.n. of Laser 2 


Laser 1 


n.i. 




n.i. 


Sum (Ss) 


none 


+ 


+ 


+ 


Diff. (Sd) 


none 


+ 


+ 





Table 1. Summary of the different possibilities of noise mea- 
surements, s.n. means shot-noise; n.i. stands for no influence. 

with the SDl and SD2 switches both in the sum position 
and changing the SD3 switch. 

3 Theoretical Model 
3.1 The phase diffusing field 

We developed a model based on Bloch equations. In the 
D2 line oi^^Rb two excited levels (out of 4 levels inside the 
Doppler broadened curve) can lead to EIT with the hy- 
perfine ground states. Inclusion of a second excited level is 
important to give a better agreement with the experimen- 
tal curves. The atom is excited by two classical fields that 
will be considered as having constant amplitudes and in- 
dependent stochastic phase fluctuations. Our model is an 
adaptation of the model of Ref.[12j for a three-level atom. 
The four levels are represented in Fig.[Tl where we made a 
correspondence of the hypothetical quantum states to the 
realistic levels of ^^Rb. The two ground states |1) and |2) 
correspond to ^^Rb 5Si/2{F = 2) and (F = 3), and the 
excited states |0') and |0} stand for the 5P3/2(-F'' = 2) and 



{F' = 3) levels, respectively. We notice that these levels 
form two A systems for atoms of two velocity classes differ- 
ing by kv ~ 64 MHz, and the EIT resonance in the room 
temperature vapor is built up from nearly equal contri- 
butions of both these A systems. Although all these four 
levels are important for a good agreement with experimen- 
tal data, for the sake of simplicity we present an outline 
of calculations for a three level system (excluding the ex- 
cited |0') level). The 3-level system can already reproduce 
many of the experimental aspects associated to the EIT 
resonance [6120) . Further inclusion of the fourth level is 
straightforward. In the following sections we will present 
numerical results for both cases. 
The laser fields are given by 



F,i{t) = Si exp [i{u!it + 4)i)]et, 



(4) 



where i = 1, 2 is a label to designate lasers 1 and 2, respec- 
tively. £i is the laser's complex amplitude, cui its frequency 
and Gj is a unit vector designating the field's polariza- 
tion. The time evolutions of the phases <j>i{t) and <j>2{t) 
are described by two independent, uncorrelated Wiener 
processes [23]. This corresponds to model the lasers as phase- 
diffusing fields, with Lorentzian lineshapes[T2]. Phase fluc- 
tuations satisfy the relations 



(5) 



where 2bj corresponds to the spectral width of laser j 
and (• ■ • ) denotes stochastic average that is taken over a 
sufficiently long time. The djk function accounts for the 
initial independence of the two lasers in our experiment, 
so they have a zero degree of correlation. 



Eo„t(0 =Ei(i)+E2(t)+z-^PW, (6) 

2ceo 
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For an optically thin sample, the output field can be write 
written as , ^ 2 f°° , , ^ , ^ 

Sqq-[UJ) =13 £q£q, I duj^q I AuJ^^, gq{uJOq)gq' {(^Oq') 

/oo 
dr exp iujT[{pq{t + r, UJOq)Pq' {t, UjQq,)) 
-00 

- {Pq{t,^nq)){Pq'{t,Uj'„q,))] + 

/OO 
drCXp iujT[{pq{t + T, UJOq)p*q: [ii^Qq')) 

where (5 is a, real constant depending on the atomic density 

-{Pq{t,LO0q)){p*q,{t,Lo'0q,))]+Cc}, (9) 

and length of the sample, and P is the complex polariza- 
tion excited in the medium given by where we can see that the power spectra are obtained 



P(f) = / da;oig(woi)pi(t,woi)exp[i(tJit-|-(/)i)] + 

— OO 

dwo25(wo2)P2(i, W02) exp [i{uj2t + 02)] ■ (7) 



from the covariance matrix of the detected intensities, 
which are then ultimately related to the covariance matrix 
for the atomic variables pi = |pi| and p2 = |p2|- We now 
present an outline for the calculation of the atomic covari- 
ance matrix for the case of a three level atom excited by 
In this expression, the inhomogeneous Doppler broaden- phase-diffusing fields. More details can be found in 

ing is given by g{ujQi), for atoms with resonance frequen- [13' especially in its Appendix B. 

cies ujQi in the laboratory reference frame. pi{t,LUoi) and 
P2(t,Ci^02) are the slowly- varying atomic coherences ex- 



3.2 Atomic polarization spectra 



cited by fields 1 and 2, respectively. The total Hamiltonian can be written as 

The detected intensities of fields 1 and 2 are given by H{t) = Hq + V{t), (10) 

Iq{t) = 2 ceo |Eojit(i) • eg(t)P, where q = 1,2. All power 

where i?o = ?i'^oi|0)(0|+?i(jJ2i|2)(2| is the free-atom Hamil- 

spectra can be obtained from the expression 

tonian and 

V{t) =-ni2iexp \i{ujit + (l)i)]\l){0\ 

/OO 
[{Iq{t + r)Iq,it))-{Iq{t + r)){Iq,m ,^^V,V/n| Hl^ 

-OO - 7iJ22 expj(w2i + 02)|2)(O|, (11) 

X exp {iuJT)dT, (8) . . . 

is the interaction Hamiltonian, with the corresponding 

Rabi frequencies for both coupling fields given by Qj. 

and we recall that the sum and difference spectra are given Since no detailed Zceman structure is considered, we took 

by Eqs. ([3]). If we discard terms of second order in (3 and both Rabi frequencies as real. We now follow straightfor- 

terms independent of r, and use Eqs. ^ to ([H]), we can ward steps to write the Bloch equations, in the Liouvillc 
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form, from (fTO|) and imposing the Rotating Wave Approx- However, products in the form (pg(t+r, wog)^^/ (i, ij^oq/)) 

imation (RWA), which results in and {pq{t + T,ujQq)pqi{t,ujQy,)) appear in equation ([9]). To 

, r-AT/ , I \^ r-AT/ , evaluate these terms it is convenient to first calculate the 

dy — exp[i]\i{Luit + (pi)\exp[i]y2{ij->2t + (p2)\Axdt 

second order correlation function 

+yodt . (12) 

HereiVi andiVa are square diagonal matrices with only ze- {G{t,t;u;oj,uJok)) = {x{t,u;oj)x^ {t,uJok)), (17) 

ros and ones, wn is a column matrix accounting for the con- , , , i i ; 

" " and then calculate 

tinuous flow of atoms through the laser beam and A(i7i, f22-, 

7^ A A X \ ■ <-u w -D1 I, J- ■ 4-1, J- • {c2it,t]UJoj,uJok)) = {G{t,t;uJoJ,uJok))- 

^,^,Al,A2,OR) IS the evolution Bloch matrix, that is j // ^ ^ j // 

function of several physical parameters: F is the total 

excited state decay rate, 7 is a decay rate for the lower 

states coherence associated to the finite interaction time, 

Aj = LOj — ujQj is the optical detuning associated to laser 

J, and 6fi = Ai — A2 is the Raman detuning. The column 

matrices containing the rapid and slowly varying elements 

of the atomic density matrix are y and x, respectively. 

They are related by the transformation 

X = exp[-iNi{uiit + (t)i)]exp[-iN2{iLJ2t + 4>2)]y- (13) 



We have special interest in the x matrix, because it con- 
tains the slowly varying atomic coherences {poi,Po2,Pi2, 
and their conjugates). To proceed with the calculations of 
the stochastic averages one must expand the exponential This and the use of the regression theorem allow one 



{xit,iJoj)}{x^t,LJok)), (18) 

where a;^ represents the hermitian conjugate of x. Finally, 
we use the regression theorem to compute {G{t + r, t; ujQj , 
<^ofc))- To obtain an equation of motion for (G(t, t\ ojqj, ujok)) 
we use the definition (fT^]) , differentiate the right-hand-side 
keeping up to second order terms in the stochastic phases 
and use resulting in 

d(G(t,t;woj, wofc)) = 
{-Ai(tJoi)(G(t,t; woj, wofc) - {G{t,t]uj()j,ujQk))A\{ujQk) 
+ 2biNi {G{t, t- LJoj ,ojok))Ni+ 2b2N2 {G{t, t; ljoj , wofc))iV2 
+ yo{x1{t,LJoj)) + {x{t,LJoj))yl}dt. (19) 



factors up to second order in the d(t)j 's and take averages to get an equation of motion for (c2 {t, t + r; uj^j , wofc)) 
using eq.®, resulting in a differential equation for (x) ^ 

d{x)^[-Ai{x)+yo]dt (14) 



— {c2(t,t + T;ujoj,uJok)) = -~Ai{c2{t,t + T;ujoj,^^ok)), 
dr 



(20) 

with which will be used in the calculation of the spectra. A 

, ■ , • AT , I, Ar2 , I, Ar2 a ^iir\ possiblc way to obtaiu a solutioii of cq. (1201) is to take its 

Ai = iNiLOi + iA'2ti'2 + OiN^ + O2JV2 — A, (15) ^ -1 H — V 

Laplace transform 

whose steady state solution is 

{x) = A-^y^. (16) 



Q{s]LOQj,ujQk) ^[s + Ai] ^{c2{t,t;uJoj,^ok)) (21) 
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4 6 8 10 12 

Analysis Frequency (MHz) 



Fig. 2. Individual laser noise for both lasers after interaction 
with the atomic sample. Circles: Laser 2 (full) and Laser 1 (hol- 
low) at high power (/ = 118 mW/cm^); Lines: Laser 2 (black) 
and Laser 1 (gray) at low power (/ = 25 mW/cm^). Resolution 
Bandwidth (RBW) 1.0 MHz, Video Bandwidth (VBW) 3 kHz. 
Each curve is an average over 100 measurements. 



where Cy(s; wq^, wofc) is the Laplace transform of (c2(t, i + 
T;ujQj,ujok)), and (c2(i, t; woj, (^ofc)) can be calculated us- 
ing the steady state solution of (fT9|l . Since the Laplace 
transform is related to the Fourier Transform by 




20 40 60 80 100 120 
Intensity ( mW/cm^ ) 




Fig. 3. Variation of the correlation coefficient with laser in- 
tensity (per beam) for an analysis frequency of 3.5 MHz. (a) 
Experiment. RBW = 1 MHz and VBW = 3 kHz. Each point 
is an average over 6 x 10'' measurements, (b) Theory (as de- 
scribed in 22} • r is the spontaneous emission decay rate from 
the excited state (_r ~ 27r x 6 MHz). Other parameters are 
Ai = A2 = 2tt X 28.6MHz, fei = 62 = O.OSr, 7 = 0.02r. 



4 Results 



'OO POO 

/(r) exp iuTdr = /(r) exp iijjTdT+ 
-00 Jq 

Jo 

^g{s^ -iuj) + g\s = iLu), (22) 



the solutions g{s; ujQj, c^ofc) will be used to obtain the final 
results of ([5]). To include the fourth level one just adds a 
new level in the Hamiltonian (eqs. (|10p ) and repeats the 
calculation. 



4.1 Experimental results 

The first step to measure the correlation between the fields 
transmitted by the atomic sample is to characterize the in- 
dividual noise spectrum of each field. Before interaction, 
the intensity noise of the each ECDL is slightly above 
the standard quantum limit (SQL). In contrast, in spite 
of their narrow linewidths (~ 1 MHz), both lasers have 
large amounts of phase noise producing a very broad back- 
ground spectrum [11 j. 
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After interaction with the Doppler broadened atomic 
sample, the transmitted fields present high degrees of in- 

^, . . . , -O.s j (a) | -o.sj (e) 

tensity fluctuations, f he intensity noise spectra extend to l o, , i.Oj 




frequencies as high as the Doppler width of the sample [9]. 



0.5- 0.5- 

O 0.0 -w. '-■--^'^^■-^^^ ..-^ . 0.0 - 



-0.5 



In the D2 hne of °'^_R&, the energy separations of all ex- i!ot , i!o 

0.5 0.5 

cited levels are smaller than the Doppler width. Thus, o 00 00 



(b) .0.5] (f) 



-0.5 



(c) -0.5^^ (g) 




Analysis Frequency (MHz) 



in a vapor cell each laser excites all the atomic transi- I'.o, i.o 
tions allowed by dipole selection rules. The line strength 
of the F — > i^'— group is proportional to the mean value 
averaged by all possible dipole transitions belonging to 
this group. Substituting in the measured values, the F — 

Fig. 4. (a)-(d) Experiment. Correlation coefficient spectra 

3 F'— group has an effective dipole moment 1.12 times ^ . „ , , . . . 2-. , ^ ■. ^ / ^ 

tor various field intensities (m mW/cm ). (a) 13, (b) 61, (c) 

greater than the F — 2 i^'— group. If both fields have r,/^ / 1\ no rm ■ / \ ■ ^^ 

° o J: 95 (^jjj iiji^ lYie gray curve in (a) gives the same measurement 

equal intensities (as is the case here), the Rabi frequency ^-^-^^^^ ^j^^ ^^p^^ ^1^^ ^-^^^^^ intensity. RBW = 

associated to Laser 2 will always be greater than the one i q mHz, VBW = 3 kHz. Each curve is an average over 100 

associated to Laser 1. As a consequence, for a sufficiently measurements, (e)-(h) Theoretical result for the correlation 

high power, the transmitted intensity fiuctuation of Laser coefficient spectra for various Rabi frequencies (J7): (e) O = 

2 wiU be greater than that of Laser 1. O.Sr, (f) n = 1.6r, (g) n = 2.0r and (h) n = 2.4r. Other 

parameters are Ai = A2 = 2tv x 28.6MHz, bi = bi = O.OSr, 



7 = 0.02r. 

The final result is a nonlinear character of the PN-to-AN 
conversion process. 



In Fig. [5] we present these noise spectra for both lasers 
at two different intensities measured after interaction with 
the atomic medium. For high power, the PN-to-AN con- 
version is considerably more efficient for the laser that is 
locked to the F = 3 ^ F'— group transition (laser 2), 

than for the laser locked to the F — 2 ^ F'— group Iran- In Fig. [3Ja) we show the experimental results for C 

sition. However, for a sufficiently low power, the efficiency as the intensity is increased, for a fixed analysis frequency 

of the PN-to-AN conversion is very low (and comparable) (w = 27r x 3.5 MHz). The first important point to notice is 

for both fields, and the noise power of Laser 2 can be lower the high degree of correlation between the two fields after 

than Laser 1 for small analysis frequencies. In this case, we the sample, which can reach absolute values above 0.6. 

have to keep in mind that absorption also plays a role, at- The other important feature is the clear transition from 

tenuating the mean field value and its fiuctuations as well, correlation to anti-correlation as the intensity is increased. 
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passing through a nearly uncorrelated situation around 
55 mW/cm^. 

We also measured the correlation spectral dependence 
which is shown in Figs. [l{a)-(d) for different intensities. 
We observe that the transition from correlation to anti- 
correlation occurs for all analysis frequencies from 2.5 MHz 
up to 14 MHz. Outside this spectral range, electronic noise 
prevents us from measuring the four power spectra nec- 
essary to evaluate C. At 14 MHz the correlation almost 
vanishes for all intensities. Nearly zero correlation is ob- 
served in the range from 50 to 60 mW/cm^, with small 
fluctuations depending on the analysis frequency. Outside 
this range, the beams are clearly correlated over all the 
observed spectrum. Without the vapor cell in the beam 
pathway (gray curve in Fig. (Ha)), the correlation goes to 
zero for all analysis frequencies, as expected for two inde- 
pendent lasers. 

The change in sign of the correlation coefficient can be 
understood as a consequence of the competition between 
two different processes occurring in a 3-level A system: 
EIT and two-photon Raman transitions (both Stokes and 
anti-Stokes). At low intensities, EIT generates intensity 
correlations between the fields, since higher intensities of 
one field lead to an increase in the transparency of the 
medium to the second one. As the intensity is increased, 
the atomic transitions are power broadened and the sys- 
tem becomes saturated, so the Raman process (where one 
photon is absorbed from one field and emitted in the 
other) dominates the atomic excitation. Since in this case 
the decrease in one field's intensity results in an increase 
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of the other's, it leads to an intensity anti-correlation be- 
tween them. This will be clarified below, when we theoret- 
ically analyze the role of the optical detuning in a 3-level 
and in a 4-level atom. 

4.2 Numerical results 

In order to understand the influence of the various physical 
parameters involved in our experimental data, we explored 
our model in different ways. First, we analyzed the situ- 
ation of a 3-level atom at rest, observing the dependence 
of correlation with the optical detuning, giving a physical 
interpretation for the change in sign of the correlation co- 
efficient. Next, we studied the effect of the Rabi frequency 
on the correlation for the same situation, and in a second 
moment we performed the Dopplcr integration, observing 
the contribution of atoms of different velocity classes. Fi- 
nally, we emphasize the contribution of the other excited 
atomic level presenting numerical results for the 4-level 
system. These results are compared to experimental data, 
showing good agreement. 

The conversion of phase- noise to amplitude-noise is the 
main source of fluctuations observed in this system, but it 
is not sufficient to explain the passage from correlation to 
anti-correlation. We can have a better understanding by 
analyzing what happens to the correlation coefficient in 
the simplified model of a 3-level system at rest. In a yl con- 
figuration, the EIT resonance occurs in a frequency win- 
dow (usually) much smaller than the natural linewidth as- 
sociated to the optical transition. In other words, the two- 
photon Raman detuning — expressed as Sr = A2 — Ai 
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Fig. 5. Numerical calculation of the correlation coefficient C 
for a 3-level atom at rest as function of analysis frequency for 
different optical detunings. Two power ranges were analyzed: 
(a)^2i = O.ir, (b)i?i — r. Other parameters are O2 = 1.12J7i, 

61 = 62 = o.osr, 7 = omr. 

following the notation of Fig.[T] — should be zero and both 
fields must be nearly resonant with the (real) atomic state 
{A2 = A\ = 0). If 8ii ^ 0, linear absorption should occur. 
On the other hand, if both fields are off resonance with 
the excited level and (J^? — 0, the Raman process prevails. 
In this case, the atom can not absorb a photon from one 
of the two fields independently from the other. Only the 
two-photon stimulated Raman process, in which a photon 
absorbed from one field is re-emitted into the other field, 
can occur with high probability. We understand that the 
competition between these two processes is the basis of 
the observed change of sign for the correlation between 
pump and probe fields. 



In Fig. Owe show results that support our arguments. 
We numerically calculated the correlation coefficient C 
as a function of analysis frequency for a 3-level atom at 
rest with constant Rabi frequency fl and different values 
of the optical detuning A. The calculations for low field 
intensity, presented in Fig. [5]Ja), give a nearly flat spec- 
trum, with a reduction in the absolute value of correlation 
for higher analysis frequencies, following the behavior ex- 
pected from the limited linewidth of the laser phase-noise. 
In these curves, we can see clearly the change from correla- 
tion to anti-correlation with an increasing detuning. This 
effect can be interpreted as the passage from the resonant 
EIT to a nonresonant Raman process, accompanied by a 
change in the photon statistics. For a higher intensity - 
Fig. [5l^b) - the anti-correlation approaches its limit even 
for higher analysis frequencies. This can be seen as a con- 
sequence of the power broadening produced by the growth 
of intensity, increasing the contribution of the Raman pro- 
cess. Therefore, two mechanisms are present. While detun- 
ing reduces the EIT process and the intensity correlation, 
power broadening increases the Raman process and the 
anti-correlation. 

A more detailed study of the effect of the field inten- 
sity can be seen in Fig. [H^a) . Here we analyze an atom at 
rest, with zero detuning. We can see the change from cor- 
relation to anticorrelation as a consequence of the increase 
in the Raman process, together with a broadening of the 
shape of the curve, demonstrated by an increase in the 
frequency for which the correlation changes sign. This can 
also be associated with power broadening of the atomic 
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Fig. 6. Numerical results of the correlation coefficient as a 
function of the analysis frequency for various Rabi frequencies 
(n) in the case of a 3-level atom at rest (a) and for a Doppler 
broadened ensemble (b). J7i — O.bF (solid), S7i = l.OF 
(dashed), f2i = 2.2r (dotted). ^2 = 1.12r2i, T f» 27r x 6 MHz 
for *^5Rb, 6i = 62 = O.OSr, 7 = 0.02r, Ai ^ A2 = 0. 

transition. Although wc can see a few similarities with 
the experimental case, such as the change of in sign of the 
correlation coefficient C with the incident intensity, the 
correlation changes rapidly with the analysis frequency, 
differently from what is observed in the experiment. 

A better agreement to experimental data is obtained in 
Fig-ISIb), where the Doppler integration was performed. 
In fact, it is well known that in the A configuration with 
co-propagating fields, atoms belonging to all different ve- 
locity classes contribute homogeneously to the signal, so 
if one is calculating the mean values, the Doppler integral 
can be avoided. However, for a phase diffusing field the 



optical detuning is a very important parameter and the 
Doppler width must be taken into account. One can see 
intuitively that a two-level atom at rest perfectly resonant 
with the field is almost insensitive to the field phase fluc- 
tuation, because it is at the maximum of the absorption 
curve. In contrast, if the atom is at the maximum slope of 
the absorption curve (or at the maximum of the dispersion 
curve), a small phase fluctuation will induce a large inten- 
sity fluctuation in the absorption profile and, as a conse- 
quence, in the transmitted light. This is the reason why 
the PN-to-AN conversion is associated to the real part of 
the atomic polarization [9II2j . Thus, the Doppler integral 
accounts for the atoms having all the possible optical de- 
tunings and the PN-to-AN conversion in an atomic vapor 
is better reproduced. Furthermore, since the only source 
of noise in the theory is the fluctuating phases of the inci- 
dent fields, this indicates that PN-to-AN conversion is the 
basic process behind our experimental observations. 

With the integration over the Doppler width, we finally 
observe a curve that has a change from correlation to anti- 
correlation, but with a profile that changes slowly with 
the analysis frequency and doesn't reach the high values 
of correlation calculated for an atom at rest. Fig. [6jb) 
still presents quantitative differences with respect to the 
experimental data. For example, we see that the passage 
to anti-correlation occurs for a broad range of analysis 
frequencies, higher than 4 MHz, while in the experiment, 
this passage occurs for smaller frequencies. As seen below, 
the inclusion of the fourth level in the model provides 
better agreement with the experimental data. 
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A better description of our system is obtained by in- 
cluding the second excited level |0') shown in Fig. [T] in the 
numerical calculation. In Fig. |3]Je)-(h) we show results for 
the correlation coefficient C as a function of the analy- 
sis frequency, for different values of Rabi frequencies, in 
the case of the 4-level system. In this situation, each tran- 
sition has a different atomic dipole moment, so the field 
intensity is parameterized by a global constant ]7, which is 
proportional to the field amplitude and the atomic dipole 
moment. We chose such a range of field intensities in or- 
der to adjust the theoretical results to the experimental 
curves. We notice that in the experimental situation the 
lasers were not locked to a real atomic transition corre- 
sponding to atoms at rest. Instead, we used a saturated 
absorption scheme to lock laser 2 on a cross-over peak, and 
the other laser was locked on the EIT resonance formed 
by the superposition of both. In this sense, the zero ve- 
locity atomic class was detuned approximately 28.6 MHz 
above the F' = 3 level. Since in the model the zero energy 
reference for the excited state is taken on the |0) state, we 
had to include an optical detuning in order to get a bet- 
ter agreement with the experimental curves. If this optical 
detuning is not considered, the anti-correlation is signifi- 
cantly reduced (in absolute value) , but the spectral feature 
of C(ti;) does not change appreciably. We can observe the 
change from correlation to anti-correlation for increasing 
Rabi frequencies, and a good qualitative agreement of the 
spectral plots, especially for higher Rabi frequencies. 

In order to compare these results with our experimen- 
tal data, we also calculated the variation of C with the 
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Rabi frequency for a fixed analysis frequency (tj = 27r x 3.6 
MHz). This is shown in Fig.[3l(b). We clearly see a change 
in sign for the correlation between pump and probe fields 
as in the experimental case, with a good agreement to the 
experimental data. 

Finally, wc now briefly comment the role of the laser 
linewidth on the sign of C{uj). In light of the previous 
analysis, the correlated fields in the EIT situation become 
anti-correlated if either the optical detuning is increased 
or the atomic transition becomes power broadened. We 
checked numerically that if the laser linewidth increases so 
much that it is comparable to (or higher than) the excited 
state decay rate F, the correlation between fields tends to 
change sign. The physical mechanism is totally analogous 
since the effect of laser broadening is to produce more side- 
bands in frequencies that are not perfectly resonant with 
the EIT transition, favoring the Raman process. We also 
confirmed that in a 3-level system this effect is more pro- 
nounced than in a 4-level atom. In other words, for the 
same laser linewidth and power, C{uj) is more negative 
in the case of a 3-level system than in the 4-level situa- 
tion. Consider that the carrier laser frequency is resonant 
with one of the excited levels. In the 3-level atom, the 
laser side bands will be far off resonance with the atomic 
transitions, producing pure Raman transitions, while in a 
4-level atom, these side bands approach the other excited 
state forming a second A system (eventually becoming res- 
onant), so both processes tend to compensate each other. 

It is important now to address some significant dif- 
ferences between the experimental and theoretical results 
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presented in Fig. ID The main discrepancy observed is that 
the theoretical curves do not show the fast decay of the 
correlation as the analysis frequency is increased. Discrep- 
ancies in the higher frequency domain can be accounted 
for by the amplifier gain and the reduction of the signal-to- 
noise ratio of our detection. Moreover, in the experimental 
curves, the Spectrum Analyzer measures the noise power 
centered at a chosen frequency and averaged over a Res- 
olution Bandwidth of 1 MHz. In the model, the analyzer 
is supposed to be ideal. 

Differences may also come from the simplicity of our 
model. In the Doppler broadened ®^Rb D2 transition, the 
excited level is composed of 4 states, namely F' = 1 — 
4, two of which (the F' = 2 and F' = 3 considered 
in the model) contribute to the A level-schemes. How- 
ever, three of the four levels contribute to the each of 
the single-photon optical transitions, which give rise to 
the PN-to-AN process. These are crucial for the Sii{uj) 
and S22{^) spectra - used for determining the correlation 
coefficient C. An important point comes from the absence 
of the F' — A level in calculating 6*22, because this is 
the strongest optical (and closed) transition and is the 
main responsible for the increase in the noise power seen 
in Fig. [5] for the Laser 2. A similar reasoning is valid for 
the F' ~ 1 level for the case of Laser 1. Furthermore, in 
Fig. [3] we see that, in the experiment, the laser intensity 
corresponding to the change in sign for the correlation co- 
efficient is much higher than the saturation intensity. In 
the theoretical case, the Rabi frequency necessary for this 
change in sign is approximately twice the saturation. We 
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remember that the model is based on the assumption of 
an optically thin medium, which is not necessarily true 
in the experiment. So, the real field intensity required to 
overcome the EIT effect and introduce anti-correlations 
via Raman transitions may be considerably higher. 



The way we model the laser noise may also give rise 
to discrepancies with the experiment. We considered that 
the lasers have perfect Lorentzian lineshapes, which have 
slow-decay spectral wings. A more realistic model for a 
diode laser, on the other hand, should include a gaus- 
sian cutoff to the Lorentzian lincshape [251 191 17] . so that 
the wings of the field spectra fall much faster than for a 
pure Lorentzian shape. This comes from the fact that the 
amplitude and phase fluctuations are strongly coupled in 
a solid state laser |26l27] . Noise spectra, from PN-to-AN 
conversion, will also depend on the shape of the phase 
noise spectra. This gaussian cutoff may account for the 
smaller absolute values of correlation obtained in the ex- 
periment, in comparison with the theoretical data, as well 
as the faster reduction of the correlation for higher anal- 
ysis frequencies. Models different from the one adopted 
here may give a better agreement to the experimental 
data, but their implementation is a much harder task 
using the present framework. Nevertheless, the simplified 
Lorentzian description already gives us a good understand- 
ing of the physical processes involved, and quite good 
agreement to the experimental data. 
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In summary, we have shown that the propagation of two 
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a broad range of analysis frequencies. We also observed the 
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transition from intensity correlation to anti-correlation as 
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